Inhomogeneous spin domain induced by the quadratic Zeeman effect in spin-1 

Bose-Einstein condensate 



Dun Zhao, 1,2 Shu- Wei Song, 1 Lin Wen, 1 Zai-Dong Li, 3 Hong-Gang Luo, 2 and Wu-Ming Liu 1 

1 Beijing National Laboratory for Condensed Matter Physics & Institute of Physics, 
Chinese Academy of Sciences, Beijing 100190, China 
2 School of Mathematics and Statistics & Center for Interdisciplinary Studies, Lanzhou 130000, China 
3 Hebei University of Technology, Tianjm 3004-01, China 

We show two kinds of inhomogeneous spin domain possessing Neel-like domain walls in spin- 
1 Bose-Einstein condensate, which are induced by the positive and negative quadratic Zeeman 
effect (QZE) respectively. In both cases, the spin density distribution is inhomogeneous and has 
zeros where the magnetization vanishes. For positive and negative QZE, the spin patterns and 
topological structures are remarkably different. Such phenomena are due to the pointwise different 
axisymmetry-breaking caused by the pointwise different population exchange between the sublevels, 
arising uniquely from the QZE. We analyze in detail the inhomogeneous domain formation and 
related experimental observations for the spin-1 8r Rb and 23 Na condensate. 



Introduction-. The experimental realization of spinor 
Bose-Einstein condensate (BEC) in an optical trap [TJ [2] 
provides a unique opportunity to study the quantum 
magnetism in a controllable way since nearly every prop- 
erty of quantum gases can be perfectly controlled in the 
laboratory [3J. For example, the spin-exchange interac- 
tion can be tuned by optical [JJ [5] or microwave [5] Fes- 
hbach resonance techniques. This is important since the 
sign of the interaction determines the ground state of 
spin-1 bosons, it is ferromagnetic (anti-ferromagnetic) if 
the interaction is negative (or positive) [7H5]. 

A more powerful way to control the spinor gases is 
to use an external magnetic field B, which leads to the 
Zeeman effects. While the linear Zeeman effect (LZE) 
is relatively trivial due to the conserved total magnetic 
quantum number, the quadratic Zeeman effect (QZE) re- 
sulted from the Zeeman energy difference in a spin-flip 
collision between different hyperfine levels[10], which is 
proportional to B 2 , becomes a key factor in determin- 
ing the properties of the spinor BEC due to its compe- 
tition with the spin-exchange interaction. The QZE not 
only leads to a number of novel ground states and spin 
structures [TUHTBj . but also affects significantly on spin 
dynamics [T9H28] . phase transitions [28H30] . symmetry 
[21 [31] [32] and the vortex states [33J , and so on. 

In spinor BEC, the QZE can be induced experimen- 
tally by laser or microwave dressing field [HI [23 CUB [33] , 
and its sign can be varied by using either positive or neg- 
ative detuning. By use of the QZE, the single-particle 
energies can be varied to control the dynamical insta- 
bilities induced by the the spin mixing collisions, these 
instabilities provide access to a rich variety of physical 
phenomena. Such controllability shows a great chance to 
manipulate spinor gas via QZE. Although many exciting 
features in spinor BEC with QZE have been investigated 
by both theory and experiment, a natural question is: in 
what extent can the QZE serve? 

In this Letter, we show that the QZE can arouse inho- 
mogeneous population exchange between the sublevels 



of spin-1 BEC, which leads to the pointwise different 
axisymmetry-breaking, and brings out the inhomogeneity 
of magnetization. We find two novel kinds of topologi- 
cally different inhomogeneous spin domain which possess 
Neel-like spin domain walls determined by unmagnetized 
areas, due to positive and negative QZE respectively. 
The spin domain can be controlled by the QZE. Further- 
more, the QZE influences deeply the spin pattern and its 
topological structure. For the positive QZE, the surface 
of space-time evolution of spin density vectors is home- 
omorphic to the quotient space obtained from two disks 
by identifying two centers. This case can in some extent 
explain the experimental observation in [17] for a spin-1 
87 Rb Bose gas. For the negative QZE, it is homeomor- 
phic to a sphere, this case is expected to be observed in 
experiment. Our results provide a profound understand- 
ing for the controllability of spinor BEC via the QZE. 

Spin-1 BEC with quadratic Zeeman effect-. We start 
from the 3-component Gross-Pitaevskii equation in a di- 
mensionless form 

ift#±i = -^<9 2 $±i + (c + ci) (|$±i| 2 + |$ | 2 ) *±i 

+(c - ci)|$ T i| 2 $±i + ci$* + (9T#±i, 

id t $o = -\d 2 x ^ + (co + d) (|$!| 2 + |$_i| 2 ) $ 

+2c 1 $*$ 1 $_ 1 +c |$o| 2 *o, (1) 

where eg and c\ are the mean-field and the spin-exchange 
interactions, p and q correspond to the LZE and QZE re- 
spectively. The energy is in units of |ci \n with n being the 
average particle density, and ci = 4-7r/i 2 (a2 — a )/(3M) 
is the spin-exchange interaction, which is ferromagnetic 
if 5i < (as 87 Rb ) and antiferromagnetic if C\ > 
(as 23 Na). Here ao and a 2 denote the s-wave scattering 
lengths of the total spin-0 and spin-2 channels respec- 
tively [7] [2]. In addition, time t, position x and $±1,0 
are in units of h/(\ci\n), h 2 / (2M\ci\n) and y/n, respec- 
tively. For 87 Rb and 23 Na, the typical particle density is 
n ~ 10 14 cm- 3 [JU1 [T3 HD], thus |ci|n is about h x 3.6 Hz 
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and h x 24.1 Hz, respectively. The time and length are 
measured in units of 10 _2 s and 10 _6 m, respectively. 

Experimentally, the LZE can be changed indepen- 
dently along the condensate length by applying a field 
gradient B' along the axis of the trapped condensate 
[TU] . taking as B' Zb where Zb is the gradient field range. 
The QZE can be changed independently by applying 
a weak external bias field Bq [TU], which is given by 
q = qB$, where q = ^ 2 /8v h f s (Sis hfs is the hyper- 
fine splitting frequency), the constant /j, is defined by 
fi = .g s /is/4/i = 700kHzG~ 2 , where g s denotes the elec- 
tron g-factor and /is is the Bohr magneton. It is found 
that q =71.65HzG~ 2 for 87 Rb and 278HzG~ 2 for 23 Na 
|10j . As mentioned above, the sign and its magnitude 
of the QZE can also be manipulated independently by 
the microwave dressing field [27l [30j [34] , which provides 
many ways to study and/or control the spin dynamics de- 
scribed by the spin density vector F — {F x ,F y , F z } [13] . 
Here F z is the longitudinal magnetization parallel to the 
external magnetic field, and one can define a transverse 
magnetization by \Fj_\ = (F 2 + F 2 ) 1 / 2 . A state with 
|-F_i_| 7^ is called a broken- axisymmetry phase ||12 . The 
polar angle determined by the interaction and the mag- 
netic field is defined by 9 = arctan(|.F_i_|/-F z ). 

In the presence of the QZE, although some simplifica- 
tion [55] and/or approximation [35- 35] have been made, 
it seems that no analytical result has been reported. Us- 
ing an ansatz with Jacobian elliptic functions, we obtain 
two kinds of analytical solutions corresponding to posi- 
tive and negative QZE respectively. 

Inhomogeneous spin domain induced by positive 
quadratic Zeeman effect-. For positive QZE, we get 
the analytical solution of Eq. (JlJ only under the fer- 
romagnetic interaction (ci < 0) with the conditions 
|ci| < c , which reads $±i = [y/—ci/(co — ci)sn(£, A) =p 
yf-q/cidnit, A)]e ix±1 , <J> = sn(£,A)e l *°, where £ = 
y/2q(x — kit), A = \f—2c\/(qS) and x±i = fa x — /i±it 
with n±i = k\/2 - q(c - 2c 1 )/(2c 1 ) - 2ci/<5 =p P, 2xo = 
Xi+X-i, 6 — (cq — Ci)/(cq+ci). The requirement |A| < 1 
gives a restriction q > —2c\/8. For 87 Rb, it becomes ap- 
proximately q > 1.982. 

In general, the above order parameter denotes a 
broken-axisymmetry phase with a periodic density dis- 
tribution, the period of the wave function can be tuned 
via q. Its spin density vector reads 

F x = \Fj_\ cospt, F y = -\F±\ sinpt, 

F z = 2 x /q/(c a -c 1 )sn^, \)dn(£, A), (2) 

where = \J — 8c 1 /(c — cf)sn 2 {^, A). The polar an- 
gle is 6 — arctan (\f~ ^f 1 dw(f'A) ) ' ^ n casc ^ ar ~ 
mor frequency of the magnetization is completely deter- 
mined by the LZE, and the polar angle is spatial and 
time-dependent, which implies the inhomogeneity of the 
pseudospin due to the QZE. 
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FIG. 1. (Color online) Inhomogeneity of the order parameter 
in spin space for the solution with q > 0. (a) and (b) show 
the cases for p — 1, q — 3, x — 1.7 with different time: t = 5.5 
for (a) and 6.5 for (b); (c) for different magnetic field (q = 9) 
in comparison to (a). The other parameters are taken from 
87 Rb atom [17] and ki = 1. 

Because an integer-spin state can be described in terms 
of the spherical harmonics Yp(s), the above order pa- 
rameter can be visualized by drawing the wave functions 
in spin space by the surface of ^(s)! 2 , as shown in Fig. 
[Tj where ^(s) = J2 m $ m F 1 m (s), and the color represents 
arg^(s) [14] , which is defined on [— 7T,7r]. Comparing 
with the ground states in a uniform system [14] . Fig. [l] 
shows the inhomogeneity of the order parameter in spin 
space. Figs, [ija) and (b) are different near the centers, 
they can not be identified via a rotation. This implies 
that the QZE can cause pointwise different axisymmetry- 
breaking in the spin evolution, which leads to the inho- 
mogeneous magnetization, as described below. Fig. [ljc) 
and (a) shows that the QZE can change the symmetry 
of the order parameter. Furthermore, the QZE can also 
change the shape near the center. 

Eq. ([2]) shows that along the propagation direction 
£ = const, such that sn(£, A) = 0, F becomes zero pe- 
riodically, which gives the unmagnetized area. When 
q -> -2ci/<5, i.e., A -> 1 (For 87 Rb, q -> 1.982, corre- 
spondingly q — > 20Hz, Bo — > 526mG), the spatial and 
time evolution of F forms a surface that looks like a 
sphere, whose south- and north-poles merge into a singu- 
lar point, corresponding to the zero point of magnetiza- 
tion. It is homeomorphic to the quotient space obtained 
from two disks by identifying two centers, as shown in 
Fig. ga). With increasing q beyond —2c\/8, two pieces 
of the spin surface separate gradually. Fig. [2^b) shows 
the time evolution of the transverse magnetization vector 
F± = {F x , F y }, where the lengths of the arrows are | F± |, 
the spin z-component density corresponds to the stripe in 
the same figure. Obviously, the magnetization is inhomo- 
geneous. The heavy yellow and red regions correspond to 
the neighbor of north pole and south pole respectively, in 
such regions, F z reaches its maximum value. The green 
and blue lines between the yellow and red regions cor- 
respond to the zeros of F, an unmagnetized area. The 
transverse ferromagnetic domain are divided by the un- 
magnetized area. While the stripes form the spin do- 
mains, the unmagnetized area forms the domain wall. 
A further analysis together with the projection vectors 
{F x , F z } and {F y , F z } confirms that the spin domain wall 
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FIG. 2. (Color online) Inhomogeneous magnetization de- 
scribed by {2}. (a) The spin density surface; (b) The spa- 
tial time-evolution of {F x ,F y } with the density plot for F z 
by color as in (a). The parameters used are p = 1, q = 2 
(for 87 Rb atom [T7 , p = 10Hz, B' z b » 0.014mG, <? = 20#2, 
Bo » 530mG), fci = 1, a: g [-5, 5],t S [0, 10]. 

has the main features of Neel wall. 

Recall the experiment in |17j . in almost quasi 1-D case, 
nearly pure spinor BECs were prepared in the unmagne- 
tized \mp — 0) phase with a large QZE. After a rapid 
quench across the quantum phase transition (of ground 
state), the transverse ferromagnetic domains form, which 
are divided by narrow unmagnetized domain walls, and 
the topological defects with unmagnetized filled cores 
have been observed. In this respect, it is constructive 
to compare the experiment with the analytical solution, 
it is easy to conclude that near q = —2ci/S, the spin den- 
sity distribution looks like dark soliton, which indicates 
exactly the topological defects with unmagnetized filled 
cores, the experimental observation is in agreement with 
the features of the solution. 

Equation ([2| also shows that the population exchange 
between the sublevels of tuf = ±1, 0, which arises only 
from the QZE, is inhomogeneous. However, it is found 
that in the propagation direction £ = the system has no 
population exchange. This unique property corresponds 
to the singular point (F = 0) of the spin density surface. 
In Fig. [3] we present pointwise magnetization as a func- 
tion of the LZE and QZE, which shows complicated and 
interesting behaviors. 

When q = —2c\/5, this solution gives a novel superpo- 
sition soliton-likc form for <I>±i and a kink soliton for <&o, 
*±i = [V-ci/(co-ci)tanh(0 T sech^)]^ 1 , 

$0 = tanh(£)e IXo . The shape of such solutions are deter- 
mined by the ratio 7 = |co/ci|. For 87 Rb, 7 ^> 1, each of 
the $±1 components behaves like a bright soliton. In this 
limit, the atom exchange between the sublevels is strong 
and thus forms two domains near the soliton propagation 
direction. On the contrary, far away from the direction, 
the atom exchange of the component m — ±1 is almost 
uniform, thus no domain can further form. As such, only 
one spin domain wall arises in this case. 

Inhomogeneous spin domain induced by negative 
quadratic Zeeman effect-. In this case, the analytical 
solution can be obtained for both ferromagnetic and anti- 
ferromagnetic interactions under the conditions |ci| < Cq. 
It is known that q > for spin-1 23 Na and 87 Rb BEC un- 
der a bias field. However, as mentioned previously, the 




FIG. 3. (Color online) Magnetization evolution in p, q of (§} 
at t = 1 for (a) x — 2; (b) x — 3.5. The parameters are the 
same as in Fig. |andpG [0,2tt], q G [1.98,10]. For 87 Rb, 
they correspond to p G [0, 63Hz] and q G [20H z, 100-ffz] or 
B'z b G [0,0-lmG] and B G [0.526G, 1.2G]. 

sign and magnitude of q can be manipulated by using 
the AC Stark shift induced by using linearly polarized 
microwaves, therefore q < is still accessible in exper- 
iments. Set v = —2q/(co — ci), Eqs. ([I]) has the so- 
lution: $±i = A ± isn(£,A )e' x±1 , $ = A n cn{£,X)e lXo 
where A$ = sj — 2A\A_\, A = \J (co + ci)v/k, and 
f = k(x - kit), x±i = k x x - fi ±1 t, 2\q = Xi + X-i, 
Mil = (kj + k 2 )/2 + 2c A 2 _ 1 -2^(c a ±c 1 )A_ 1 ±civTP 
for given k, ki, and the coefficients satisfy A1+A-1 = %/v 
and A\A-\ < 0. We have the restriction > q > 
— A^ x {cq — ci)/2. In this solution, while q = gives 
a periodic polar state, q = — A 2 _ 1 (c — Ci)/2 (if Ai = 0) 
presents a ferromagnetic phase. For a median q, this so- 
lution indicates a broken-axisymmetry phase. This sug- 
gests that by tuning q, it is possible to get polar phase 
or ferromagnetic phase from the broken-axisymmetry 
phase. In this case, the population transfer between dif- 
ferent components is determined by the QZE, tuning q 
will lead to the population transfer. In particular, for 
k = yj (co + c\)v (A = 1), the populations of all three 
components can be obtained analytically [41 , from which 
two kinds of population transfer are identified: one is the 
population exchange between sublevel m = ±1 and their 
backgrounds, the other is the exchange between sublevels 
m = and m = ±1. 

When both A±\ 7^ 0, the solution represents a broken- 
axisymmetry phase with spin density vector given by 

F x = ~2yJ~A 1 A_ 1 vsn{^, X)cn(£, A) cos Qt, 

F y = 2\J -A\A-\Vsn{ii, A)cn(^, A) sin tit, 

F z = (is-2VvA_ 1 ) S n 2 (£,\), (3) 

where f2 = (p — C\V + 2c\A_\y/v) gives the Larmor fre- 
quency. The polar angle 6 = arctan( ^= _^~ ' ) ■ 
Same as the solution for q > 0, along the propagation 
direction £ = 0, the spin density vector becomes zero 
periodically, which gives unmagnetized regime, and the 
magnetization is inhomogeneous. 

Taking 23 Na as an example, Fig. |4^a) gives the spin 
density surface for x G [—5,5], t G [0,10] in which the 
south pole of the ellipsoid is the zero of F at which the 
magnetization vanishes. In Fig. [4^b), the stripes form 
the spin domains and the transverse ferromagnetic do- 
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FIG. 4. (Color online) Inhomogeneous magnetization de- 
scribed by Eq.p}. (a) The spin density surface; (b) The time- 
evolution of the transverse magnetization vector {F x ,F y } 
with the density plot for F z by color as in (a). The parameter 
used are p = 1, q = -0.47 (for 23 Na atom 30!, p = 25Hz, 
B'z b « 0.035mG, q = -ll.&Hz), ki = 1, A-i = -0.2, k = 2. 

mains. A similar analysis as for q > 0, the formed spin 
domain wall is also a Neel-like wall. 

In Fig. ^b), while the red regions correspond to points 
near the north pole, the yellow lines correspond to the 
south pole where spin densities of all three components 
trend to zero. Other than the case of q > 0, a spin 
domain is divided into two transverse ferromagnetic do- 
mains, with the domain wall being the center line of the 
red region, where = but \F Z \ reaches its maxi- 

mum value. In the two sides of the domain walls, the 
magnetization direction is opposite. However, when q = 
— fc 2 <5/2, the solution becomes $±i = A±\ tanh£e IX±1 , 
$o = A sech £e IXo . The spin density surface still looks 
like an ellipsoid, but |Fj_| behaves like a split bright 1- 
soliton, there are two separated transverse ferromagnetic 
domains only concentrating around two lines correspond- 
ing to the splitting peaks of the soliton. The topological 
structure of the spin configuration in p, q given by Eq. 
([3]) is different from that of Eq. ([2]), as shown in Fig. [5] 

In the recent experiment |30j . by rapidly switched the 
QZE from positive to negative values, the dynamics of an 
antifcrromagnetic sodium BEC quenched across a quan- 
tum phase transition has been observed. Here we point 
out that, if preparing pure spinor BECs in the \mp = 0) 
phase at a large negative QZE, by the same technique as 
in |30j . after a rapid quench but not going across 9 = 0, 
the inhomogeneous magnetization given by the solution 
could be observed by probing the Larmor precession |42j . 

Conclusion-. Based on analytical solutions, we have 
displayed two kinds of inhomogeneous spin domain pos- 
sessing Neel-like domain walls in spin-1 BEC in the 
presence of positive and negative QZE respectively. 
Some novel inhomogeneous phenomena, which are due 
to the pointwise different axisymmetry-breaking arising 
uniquely from the QZE, have been shown. We present 
detailed description for the inhomogeneous spin pattern 
and the topological structures. The transverse magneti- 
zation with positive QZE is found to be consistent with 
the experimental observation. The inhomogeneous mag- 
netization for negative QZE is expected to be observed 
in experiment using the method in [30 , 42 . These results 




FIG. 5. (Color online) Magnetization evolution in p, q of 
Eq.Q for different positions at t = 1. (a) x = 2; (b) 
x = 3.3. Here p G [0, 10], q G [—1.88,0], which correspond 
tope [0,247 Hz] (B'zb « 0.35mG ) and q G [-46.5#*, 0]. 
The other parameters are same as in Fig. [4] 

display some novel phenomena in quantum magnetism, 
and provide a new prospect for manipulation of spinor 
BEC via the QZE. Our discussion is possible to be ex- 
tended to spin-2 BEC case. 
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